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ON POINT TRANSFORMATIONS OF LINEAR 
EQUATIONS OF MAXIMAL SYMMETRY 


JC NDOGMO 

Abstract. An effective method for generating linear equations of 
maximal symmetry in their much general normal form is obtained. 
In the said normal form, the coefficients of the equation are differ¬ 
ential functions of the coefficient of the term of third highest order. 
As a result, an explicit expression for the point transformation re¬ 
ducing the equation to its canonical form is obtained, and a simple 
formula for the expression of the general solution in terms of those 
of the second-order source equation is recovered. New expressions 
for the general solution are also obtained, as well as a direct proof 
of the fact that a linear equation is iterative if and only if it is 
reducible to the canonical form by a point transformation. New 
classes of solvable equations parameterized by arbitrary functions 
are also derived, together with simple algebraic expressions for the 
corresponding general solution. 


1. Introduction 

In a short paper published by Krause and Michel [1] in 1988 certain 
specihc properties of linear equations of maximal symmetry were estab¬ 
lished. In particular, that paper shows that such equations are precisely 
the iterative ones, and equivalently those which can be reduced to the 
canonical form = 0 by a point transformation. However that short 
paper left a number of important questions unanswered. It does not 
provide for instance any expression for the point transformation map¬ 
ping a given equation of maximal symmetry to the canonical form. On 
the other hand it also makes use without proof of a number of impor¬ 
tant results which were visibly already known but which are difficult 
to hnd in the literature. 

Almost at the same time the problem of generation of equations of 
maximal symmetry was considered by Mahomed and Leach [2] who 
found an algorithm for obtaining expressions for the most general nor¬ 
mal form of this equations based on the direct computation of the 
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symmetry algebra. However, computations with this algorithm remain 
quite tedious and the authors managed to provide a general expression 
for the linear equations of maximal symmetry only up to the order 
eight. 

A more direct algorithm for generating this class of equations based 
on the simple fact that they are iterative was proposed recently in 
[3]. Nevertheless, some of the main results of that paper concerning 
in particular the generation and the point transformations of the class 
of linear equations of maximal symmetry still have room for improve¬ 
ments. 

In this paper, we provide a much simpler differential operator than 
that provided in [3] for generating linear iterative equations of a general 
order. As a result, we give a more direct proof than that of Krause and 
Michel [1] to the fact that linear equations reducible by an invertible 
point transformation to the form = 0, and which will henceforth be 
refereed to as the canonical form, are precisely those which are iterative. 
We also establish several results concerning the solutions of this class 
of equations, and in particular their transformation to canonical form. 
In contrast to the very well-known paper by Ermakov |1] who managed 
to hnd only some very specihc cases of a very restricted class for which 
the second-order source equation is solvable, we provide full classes 
parameterized by arbitrary functions and for which the solution of the 
second-order source equation, and consequently of the entire class of 
equation of maximal symmetry, is given by simple formulas. 

2. Iterations of linear equations 

Let r 7^ 0 and s be two smooth functions of x, and consider the dif¬ 
ferential operator d' = -|- s. We shall often denote by F[xi ,..., Xm\ 

a differential function of the variable xi,..., Xm- Linear iterative equa¬ 
tions are the iterations = 0 of the hrst-order linear ordinary 

differential equation \[^[?/] = ry' + sy = 0, given by 

^n^y] ^ for n > 1 with = I, (2.1) 

where I is the identity operator. A linear iterative equation of a 
general order n thus has the form 

= a:; !/<”-*>+!/'"-"*+••■+AT'j/'+if; J/ = o. (2.2) 

Setting 

= 0, for j < 0 or j > m, and = 1, for m = j = 0 (2.3) 
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and applying (12.ip shows that the coefficients in the general expres¬ 
sion fl2.2p of an iterative eqnation satisfy the recnrrence relation 

Ki^ = rKi_^ + for 0 < j < n, and n > 1. (2.4) 

Moreover, setting j = 0 or j = n in (12. 4 p shows by indnction on n that 

for all n > 1, (2.5) 

and applying fl2.4p recnrsively and using the conventions set in fl2.3p 
give the new recurrence relation 

n 

KI = ^ kIzI, ioY j = 0,... ,n and n > 1. (2.6) 

k=j 


We note that fl2.6p provides an algorithm for the computation of the 
coefficients in terms of the parameters r and s of the source equation 
and the operator 'h, and the resulting formula has effectively been 
obtained in [3l Theorem 2.2]. Moreover, it is of course also possible to 
compute the directly in terms of n and the parameters r and s, and 
for j = 1, 2 we hnd that 




ns + \ r' 

' 2 


3n-5 


2 pW + lg) (3.r' + rr"+ ^ 


.'2 


(2.7a) 

(2.7b) 


If we divide through the general nth order linear iterative equation 
T”)?/] = 0 in fl2.2p by 77° = r"', it takes the form 


0 = ?/(") + Bl + • ■ ■ + + --- + B^y, (2.8a) 

Bi^ = KHr-. (2.8b) 


It is clear that this equation represents the standard form of the general 
linear iterative equation with leading coefficient one. Moreover, the well 
known change of the dependent variable y i—)■ y exp Bl^{t) dt^ 

maps (I2.8ap into its normal form in which the coefficient of the term 
of second highest order has vanished. This transformation however 
simply amounts to setting = 0, i.e. = 0. Therefore, for given 
parameters r and s of the operator T, an nth order linear equation in 
normal form 


y(n) + ^2 „(»-2) + ... + Ai + . . . !, = 0 


(2.9a) 
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is iterative if and only if 


(2<j<n), (2.9b) 

where is given by (I2.6p . It follows from (I2.7ap that the reqnirement 
that = 0 holds is eqnivalent to having 

s = -^(n-l)r', (2.10) 

and this shows in particnlar that any iterative eqnation in normal form 
can be expressed in terms of the parameter r alone, i.e. depends on 
a single arbitrary fnnction. Clearly, the coefficients can also be 
expressed solely in terms of n, r and the derivatives of r. For instance, 
by setting for any given fnnction 


Al = 




[e'(x)]2 - 2e(x)e"(x) 


4[^(x)]2 

it follows from (]2.7bp that in (12. 9 p we have 


( 2 . 11 ) 


= A{r)^ and more generally 


( 2 . 12 ) 


In fact, as already noted in [2], [3], the coefficients depend only on 
the fnnction A{r) = A 2 and its derivatives. For simplicity, it will often 
be convenient to denote the coefficient A 2 of the term of third highest 
order in fl2.9ap simply by q. 

Having noted that the coefficients of every linear eqnation of maxi¬ 
mal symmetry in normal form depend solely on q and its derivatives, 
an important problem considered in [2] was to hnd a linear ordinary 
differential operator F„[|/] depending solely on q and its derivatives, and 
which generates the most general form of the linear nth order eqnation 
of maximal symmetry in the dependent variable y = y{x). In a recent 
paper |3], it was established that for an arbitrary parameter r of the 
sonrce eqnation, the operator 


$ 


n 



Kk=0 


(2.13) 


generates the linear iterative eqnation of an arbitrary order n in nor¬ 
mal form and in its most general form fl2.9p . Therefore, althongh the 
operator fl2.13p and the eqnation it generates depend explicitly on r 
(and not on q = H 2 ) and its derivatives, on the basis of the resnlt of [1] 
identifying linear iterative eqnations with linear equations of maximal 
symmetry, the two equations F„[|/] = 0 and ^n[y] = 0 should always 
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be the same for all n, although such an operator r„[?/] has not yet been 
found. This is due in part to the fact that a general expression for 
Ai^ = A^[q] is not available for all n and j. Nonetheless, it is naturally 
possible to make use of the differential operator to generate directly 
a linear equation of maximal symmetry of the general form fl2.9ap in 
which the coefficients depend only on q and its derivatives. 

Indeed, for any value of n > 2 it follows from the equations fl2.6p and 
fi 2 . 9 p - fi 2 .r 2 p and a simple induction on j > 2 that each coefficient in 
fl2.9p depends linearly on Moreover, it follows from fl2.lip - fl2.12p 
that 

= ]j3-2 ;= Fj[r, q], for j > 2, (2.14) 


where = d/dx. Therefore, applying the substitution fl2.14p to ^n[y\ 
yields the desired equation, that is the linear equation of the form 
fl2.9p . in which = v4^[q] depends only on q and its derivatives. More 
formally, the resulting differential operator can be represented as 


4> 


r 

n 



iCi=0, 


r A) = Fj [r, q] 

b >2 


= $ 


n 


= Fj[r,q] 

b >2 


(2.15) 


For example, for n = 3 or 4, evaluating $„[?/] yields the following 
expressions directly in terms of r and its derivatives alone. 


y (r'^ — 2rrV"-|-r^r®) y'(A — 2rr") /o-, 

Hy] = - 5 - - + ——^^ 


(2.16a) 


44|!/] = 


3y ( 27 A — QSrAr" 24r^rV^^^ -|- 4r^ (7^'“^ — 


16 r^ 

5y' (r'^ — 2rr'r" + ^ 5 {A — 2rr") y" ^ ^ 4 ^ 


2r2 


(2.16b) 


However, if in addition we also apply to these expressions for ^n[y] the 
substitution fl2.14p . which amounts to applying directly the operator 
•FO to y we obtain 

4)^[y]=2q'y + 4qy' + y'" (2.17a) 

^l[y] =3y (3q2 + q") + lOy'q' + lOq/ + (2.17b) 

Comparing this with the known expressions for linear equations of max¬ 
imal symmetry expressed solely in terms of q and its derivatives [21 [3] 
shows that 4)3 and and 4>4 yield indeed the indicated expressions. 

If we set r = V? for a certain nonzero function u, the expression for 
A{r) = A{v?) in fl 2 . 1 ip is much simpler and reduces to —u''/u. Setting 
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q = A{r) is thus equivalent to letting u be a solution of the equation 


y" + qy = 0, 


(2.18) 


which is referred to as the second-order source equation for (I2.9|l . 
Thus, if we express ^n[y\ with r replaced by the substitution rule 
for udl similar to that given for in fl2.14p takes the much simpler 
form 


„(i) = ]ji-2 ^ > 2. (2.19) 


Denote by 0“ the operator in which r is replaced with and to 
which the substitution (I2.19h is applied. That is, 


0 “ = $. 


r = u 

^ Hj[u,qi\, j>2 


( 2 . 20 ) 


That is, 0“ is just <1>„ as in fl2.13l) . in which the source parameter r 
has been replaced by and to which the substitution fl2.19p is then 
applied. Then for the same reasons that <1)0[|/, q] generates the most 
general form of linear iterative equations, 0“[y, q] also does the same. 
However, in view of the most simpler substitution rule fl2.19p . the com¬ 
putational cost for generating these equations is much lower using 0“ 
rather than $(). Indeed for arbitrary values of q, generating equations 
of maximal symmetry of order greater than ten has been up to now 
a very tedious task using the original algorithm of [2], but it is now 
easy to generate such equations up to the order thirty in most stan¬ 
dard computers using computing systems such as mathematica. For 
instance, the coefficient A\^ of y which is the largest expression in size 
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in the equation Qi^[y\ = 0 is given by 

=14(2(52022476800q®q' + 2132810240q'^ + 6656237568q^q(^) 

+ 341232100q'3q(^) + 1024q^ (73853676q'q" + 286397q(®)) 

+ 32q3(2646561024q'3 + 207959056q"q(3) + 123346720q'q(^) 

+ 184297q(^)) + 2q'2 (863911980q"q(^) + 784597q(^)) 

+ 8q2(2094143648q'^q(^) + 33615550q®q(^) + 22652990q"q(®) 

+ q' (2835404512q"2 + 9750858q(®)) + 7323q(^)) 

+ 13(457296[q(^)]^ + q^^^ (740970q"^ + 2717q^®)) + 2068q^^)q(^) 

+ 2q(3) (I100540q"q(^) + 591q(®)) + 491q"q(^)) + q'(780095196q"^ 
+ 8478360[q(^)]2 + 14231306q(^)q(^) + 8279146q"q(®) + 1793q(^°)) 
+ 2q(l2696730880q'^q" + 151883460q'^q(®) + 2519682q(^)q(®) 

+ 2q(3) (257461378q"2 + 913529q(®)) + 937436q"q(^) + 139q(^^) 

+ 4q'(95070193[q(3)]^ + 152580285q"q(^) + 79587q(®)))) + 

Using the operator 0“, we shall give in the next section a much sim¬ 
pler and direct proof than that of [T] that a linear equation is iterative 
if and only if it is reducible by a point transformation to the canonical 
form. To close this section, we note that the Wronskian of any two 
linearly independent solutions u and v of the source equation fl2.18p is 
a nonzero constant and will be normalized to one. Nontrivial solutions 
of the second-order equation fl2.18p are not known for arbitrary values 
of the coefficient q. 


3. Point transformations 

The equivalence group of linear nth order equations in normal form 
fl2.9ap is well known to be given by invertible point transformations of 
the form 

x = f{z), y = X[f{x)]—w, (3.1) 

where / is an arbitrary locally invertible function and A an arbitrary 
nonzero constant [3 El [3- In other words, an equation of the form 
fl2.9ap is reducible to the canonical form if and only if there exist a 
transformation of the form fl3.ip that maps such an equation to the 
canonical form. 

Let us denote by 

(-3;« + 2e'^(3)) 

25'2 


mu ) 


(3.2) 













the Schwarzian derivative of the function ^ = ^(^). By studying the 
expression of the source parameter of the transformed equation un¬ 
der equivalence transformations, a simple characterization of the point 
transformation that maps an iterative equation to the canonical form 
was found in [Sj Theorem4.3]. This result states that a point trans¬ 
formation reduces a given iterative equation, which without loss of 
generality may be assumed to be of the form fl2.9ap . to the canonical 
form w^'^\z) = 0 if and only if it is of the form fl3.ip . where / is the 
inverse of the function z = h{x) satisfying 

Al(x) = -S{h){x), or equivalently, ^(r)(a;) = -S{h){x), (3.3) 

on account of fl 2 . 12 p . assuming that r is the source parameter of the 
equation. By making use of this result of [3], we can now provide a 
more direct proof than that given in [1] for the following result. 

Theorem 1. A linear equation is iterative if and only if it can he 
reduced to the canonical form by an invertible point transformation. 

Proof. As usual one may assume that the iterative equation is in its nor¬ 
mal form and has source parameter r. On the basis of the above stated 
result from Theorem 4.3 of [3], it follows that the transformation fl3.1|) 
where / is the inverse of the function h = f — maps the iterative equa¬ 
tion to its canonical form, as the latter expression for h solves (13.3p . 
More directly, one can prove that the class of iterative equations and 
that of equations reducible by an invertible point transformation to the 
canonical form are the same. Indeed, let fln[f, A, x, y] = fln[f, A] repre¬ 
sent the element of the equivalence group of linear nth order equations 
A = An[y] = 0 in normal form, acting on the space of independent 
variable x and dependent variable ?/, and given by fl3.ip . Denoting 
by fln[f, A] ■ A = 0 the transformed equation, it follows that we have 
exactly 

$„[|/] = Kn{X,r) Qri 

where 

Kn{X,r) = 1/ (A[r(a:)](”+^)/^) . (3.4b) 

In other words, generating a linear iterative equation in y = y{x) for 
any given value of r using the differential operator and transforming 
the canonical equation wA'>(^z) = 0 in the new variables x and y us¬ 
ing the point transformation operator Kn{X,r) [/ A, 2 :, tc] yield 

exactly the same equations, and this proves the result. □ 


dx 

—,A, z,w 
r 


■w^'^\z) 


(3.4a) 
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In [T] Krause and Michel proved this theorem indirectly as a con¬ 
sequence of equivalence relations, by proving the equivalence between 
equations of maximal symmetry and iterative equations on one hand, 
and between equations reducible by invertible point transformations 
and equations of maximal symmetry on the other hand. The part in 
that proof stating that every equation reducible to the canonical form 
has maximal symmetry is due to Lie [8]. The direct proof given here 
will allow us to construct various point transformations of practical im¬ 
portance for iterative equations as well as various forms of their general 
solutions. 

It should hrst be noted that the right hand side of fl3.4ap gives an¬ 
other methods for generating linear iterative equations with source pa¬ 
rameter r. In practice however, linear iterative equations arise generally 
with no reference to any source parameter but are expressed solely in 
terms of the coefficient q = ^ 42 ( 0 ;) and its derivatives. In this case, on 
the basis of a remark made in the previous section, a solution h to the 
equation q = ^S{h){x) in fl3.3p is simply given h = J ^ where u is a 
solution of the second-order source equation fl2.18p . Consequently, if 
we denote by the transformation operator given by 


n 


u 

n 


dx 


Kn{X,U^) fin 


dx - 


b '>2 


(3.5) 

where Hj[uc\\ = Di ^(—uq) as in 02.191) . then 0“[2/] and • w^'^'>{z) 
generate exactly the same linear iterative equation expressed solely in 
terms of q and its derivatives. 

On the other hand, since the transformation operator [J A, z, w'j 
maps the reduced equation to the most general form 02.9ap of the iter¬ 
ative equation, its explicit expression can be used to derive the general 
solution of 02.9ap . Indeed, under this operator we have 


y = ^'1^, 

^ A 



(3.6) 


Consequently, n linearly independent solutions to the general nth order 
iterative equation fl2.9aD are given by 

A; = 0,...,n-1. (3.7) 



In particular, if u and v are two linearly independent solutions of the 
source equation fl2.18p . then / ^ = b Consequently, in terms of u and 
V, the n linearly independent solutions to (I2.9ap in (13.7p above can be 

























10 


rewritten as 

Uk = , A; = 0,..., n — 1. (3.8) 

Formula 03.81) is well known and was cited without proof in [T], and it 
is an important result for which we have not been able to hnd the proof 
in the recent literature. Moreover, formula 03.7p established above 
provides a much stronger result by showing that linearly independent 
solutions of 02.9al) can be expressed solely in terms of a single nonzero 
solution of the second-order source equation. Indeed, the solutions yk 
in 03.71) are clearly linearly independent as their Wronskian equals the 
nonzero constant 11^=1 7•• 

Theorem 2. A linear equation is iterative if and only if it has n linearly 
independent solutions yk of the form 03.81) . where u and v are two 
linearly independent solutions of the corresponding second-order source 
equation 02.18p . 

Proof. The fact that a linear iterative equation has linearly independent 
solutions of the stated form is established in 03.81) . Conversely, if a 
linear equation has n linearly independent solutions of the stated form 
03.81) . then since the second-order source equation is reducible to the 
canonical form = 0 by a point transformation, without loss of 
generality we may assume that such a transformation reduces m to 1 and 
V to Consequently, the corresponding solutions of the linear equation 
are polynomials of degree at most n — 1, and thus the transformed 
equation is in canonical form. It then clearly follows from Theorem [1] 
that the equation is iterative. □ 

4. Concluding remarks 

The results obtained in this paper show amongst others that linear 
equations of maximal symmetry are highly solvable, because that is 
indeed the case for the second-order source equation fl2.18l) whose so¬ 
lutions completely determine those of the linear equations of maximal 
symmetry. In a very popular paper [1] published in Russian in 1880 and 
recently translated into English [9], Ermakov stated that the majority 
of second-order linear homogeneous ordinary differential equations for 
which it is possible to find conditions for their solvability are of the 
form 

-f a 2 X + af)y" + {a 4 ^x + a^)y' a^y = 0, (4.1) 

where the aj are some arbitrary constants. He then moved on in the 
paper to obtain some very specihc cases of equations which are solvable 
from this class, together with their general solutions. A couple of years 
later in the same decade Hill [10] considered in the study of lunar 
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stability the most general form of the second-order linear equation in 
normal form, i.e. in the form fl2.18p . but in which the coefficient q is a 
periodic function. Hill’s equation as well as its important variants such 
as Meissner Equation and Mathieu Equation have been been amply 
studied m, and it is well known that its solutions and their related 
properties can be described by means of Floquet theory. This solutions 
can also be expressed in terms of Hill’s determinant |T2]. 

The results of this paper extend even more the scope of solvability of 
linear second-order equations obtained by Ermakov. Indeed, it follows 
from fl3.7p and the condition r = v? relating the parameter of the source 
equation and the solution of the corresponding second-order equation, 
that for arbitrary values r = r{x) of the source parameter, two linearly 
independent solutions of the second-order equation of the form 

y'' + A{r)y = 0 (4.2) 

are given by yk = \A’(/ for fc = 0,1. In passing, we note that 
n linearly independent solutions of the corresponding nth iteration of 
(14.21) can also be obtained from (13.71) through the substitution r = u^. 
To make a much direct comparison with the class of equations (14.11) 
considered by Ermakov, we note that linearly independent solutions 
can also be found for the standard form of (14.21) which in terms of the 
arbitrary coefficient B = B{x) of y' takes the form 

y" + By' + ^(4M(r) + B^ + 2B')y = 0. (4.3) 

The two linearly independent solutions of (14. 4 p are given by 

» = J=0.1. (4.4) 

We note that as opposed to the very specihc cases of solvable equations 
found by Ermakov from the much restricted class of equations (14.ip . 
not only (14. 4 p depends roughly speaking on two arbitrary arbitrary 
functions, but also their solution is given by a very simple formula. 
One difficulty with the type of solutions found here for equations of 
maximal symmetry is however that in practice their coefficients are 
not generally given explicitly in terms if the source parameter r, but 
solely in terms of q and its derivatives, and solving A{r) = q for r is 
equivalent to solving the corresponding source equation (I2.18P through 
the linearizing transformation r = v?. 
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